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Chapter 1 

Relativistic Fluid Dynamics 



In this chapter I review the basic theory of fluid kinematics and dynamics (without dis- 
sipation) in relativistic spacetime. The classic paper in this field is Ellis' 1971 review |jT|]. 
That paper is at a more advanced level than these lectures. For a basic introduction to 
tensors, relativity and fluids, see for example [0. 

I use units in which the speed of light in vacuum, Einstein's gravitational constant 
and Boltzmann's constant are all one: 

c = SttG = k = l 

I use A = B to denote equality of A and B in an instantaneous orthonormal frame at a 
point (defined below). 

1.1 Brief Review of Relativity 

The observed universe is a 4 dimensional spacetime. Physical laws should be expressible as 
equations in spacetime that are independent of the observer. Together with experimental 
and observational evidence, and further principles, this leads to Einstein's relativity theory 
- special relativity in the case where the gravitational field may be neglected, and general 
relativity when gravity is incorporated. 

Local coordinates, which are typically based on observers, are usually chosen so that 
x° is a time parameter and are space coordinates. A change of coordinates (or of 
observers) is 



X = [X ,x 



°-'^-(t,f) ^ = (x°',xO = (t',f' ) (1.1) 



Physical laws should then be invariant under such transformations. This means that these 
laws are expressible in terms of tensor fields and tensor-derivatives. Tensors have different 
types (r, s), but they all transform linearly under ( |1.1|) . The simplest example is a scalar, 
which is invariant. Using the chain rule, the tranformation of the coordinate differentials 
is seen to be linear: 

dx" = — — dx" = — — dx°' 
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Extending this to partial derivatives of scalars and generalising, we are led to the trans- 
formation properties of tensors in general: 



(0, 0) scalar 


/ - 


f 




(1,0) vector 


- 






(0, 1) covector 




. , - 

" ax"' " 




(1, 1) tensor 


rjia 

J- p 


^' dx» 





(r,s) tensor J"^ ' "'-/3,.../3, 

It follows that if a tensor vanishes in one coordinate frame, it vanishes in all frames. 
Consequently, if two tensors are equal in one frame, they are equal in all frames. 

Fields and equations that transform according to ( |1.2D are called tensorial or covariant. 
Restricted covariance arises when the class of allowable coordinate systems is restricted. 
In special relativity (flat spacetime), one can choose orthonormal coordinates which 
correspond to inertial observers, and if is required to be also orthonormal, then 

^ = A"'„ ^ x"' = A"' + C (1.3) 

where A, C are constants and A is a Lorentz matrix. In other words, special relativity says 
that the laws of physics (leaving aside gravity) are invariant under Lorentz transformations 
that connect any inertial observers in relative motion. Under this restriction, the partial 
derivatives of tensors transform according to ( |1.2| ), i.e. they are Lorentz covariant. We 
use the notation 

J-...,,,^a,J-...,^ Aj-..., (1.4) 

for partial derivatives. Thus in special relativity, physical laws are expressed in orthonor- 
mal coordinates as PDE's; for example the Klein-Gordon equation for a massless scalar 
field is 

m = r/"^a„(9^^ = (1.5) 

where 

77„^ = diag(-l,l,l,l) = r/"^ (1.6) 

are the orthonormal components of the metric tensor. 

The metric gap of any (in general curved) spacetime determines the spacetime interval 
between events, the scalar product of vectors, and the raising and lowering of indices on 
general tensors: 

ds'^ = gafsdx'^dx'^ (1.7) 

U-V = QapU'^v'^ = u"Va = UaV" (1.8) 

^",3m = 9'"'9p.Ju% , etc. (1.9) 
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where the inverse metric is defined by g'^^g^p = 6°'i3 ■ 

The metric is a symmetric tensor. For any rank-2 tensor, we can define its covariant 
symmetric and skew parts: 

= I iy^P + Vp^) , V[„;3] = \ iy^p - Vp^) (1.10) 

so that gap = g{al3)- 

At any point (or event) P, an observer can choose coordinates that bring gapiP) 
into orthonormal form. I will call such a coordinate system an instantaneous orthonormal 
frame (lOF), characterised by 

gal3 = Vo/3 ^ gaf^iP^^^ = r]al3 (l-H) 

At each event along the observer's worldline, the lOF is in general different. In fact an 
lOF is orthonormal in a neighbourhood of the original point P if and only if the spacetime 
is locally fiat. 

In curved spacetime, the partial derivative ( [L.3D is not covariant (except when J is a 
scalar). The metric defines a connection that 'corrects' for the variations in the coordinate 
basis (equivalently, that provides a rule for parallel transport of vectors): 

r"/3<T = ^g"^^ {gfii3,a + ga^l,f3 - gfia,ti) = r"(/3<x) (1-12) 

The connection, which is not a tensor since it corrects for non-tensorial variations, defines 
a covariant derivative 

f;a f,a 

J" ...j3-a = J"' ■■■/3,cr + r"^o-^^ ■■■/3 + ■•• — •■■ — r'^/3o-J° ...^ (1-13) 

We also write Vo-J"".../? for the covariant derivative. One can always find an lOF at any 
event P such that the connection vanishes at P: 



= ^ J"-.../3;^ = (1.14) 

From now on, any lOF will be assumed to have this property. 

The connection also defines a covariant measure of spacetime curvature - the Riemann 
tensor: 

DO! -na I pa i pa per pa per / 1 i r\ 

Curvature is fundamentally refiected in the non-commutation of covariant derivative^, 
as given by the Ricci identity 

"'^ except for scalars: /;[a/3] — 0. 



4 



and its generalisations for higher rank tensors. The trace-free part of the Riemann tensor 
is the Weyl tensor C^'i^^y, which represents the 'free' gravitational field and describes 
gravity waves, while the trace gives the Ricci tensor and Ricci scalar 

which are determined by the mass-energy-momentum distribution via Einstein's field 
equations 

Rap — 2^9aP — Taf3 (1-18) 

where T^p is the energy-momentum tensor, discussed below. The Ricci tensor obeys the 
contracted Bianchi identity 

=0 (1.19) 



1.2 Fluid Kinematics 

Consider the motion of a particle with rest mass m. An observer records the particle's 
history - its worldhne - as x° = {t,x^{t)). We need a covariant (observer-independent) 
description of the worldline and velocity of the particle. If m > 0, then along the worldline 
ds^ < (the particle moves slower than light). If r is the time recorded by a clock 
comoving with the particle, the worldline is given by = x°(r), independently of any 
observer. The covariant comoving time is called the proper time. In an lOF ds^ = —dr"^. 
Since both sides of this equation are tensors (scalars), the equation holds in any frame, 
and at all points along the worldline, i.e. ds'^ — —dr"^. The kinematics of the particle are 
covariantly described by the 4-velocity 

dx" 



m" = — u''Ua = -l (1.20) 



and the 4-acceleration 

It" = U",f3u'^ (1.21) 

where u°'Ua = 0. The particle moves in free-fall, subject to no non-gravitational forces, 
if and only if iia = 0, in which case its worldline is a (timelike) geodesic. In the observer's 
lOF 

= o'(t.)(l, f ) = l{l,v) , 7(t') = (1 - vY'^' = § (1-22) 

where t is the observer's proper time at that point, and v is the measured velocity of the 
particle. 

If m = 0, the particle (photon or massless neutrino or graviton) moves at the speed 
of light, and along its worldline ds^ = 0, so that proper time cannot parametrise the 
worldline. In the lOF of an observer m", the light ray has angular frequency u and wave 
vector k (where \k\ = lu), with phase (j) = k ■ x — ut, so that 

0,a = (-<^,^) and (f)^a(f)'°' ^ 

Now the phase is a covariant scalar, and its gradient is a covariant null vector, which we 
call the 4-wave vector, and which is geodesic: 

ka^(f),a and kak"" ^ k'^.pk'^^0 (1.23) 
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From the above, in the observer's lOF, uj = —kaU'^ = 0. This gives a covariant expression 
for the redshift between events E ('emitter') and R ('receiver') along a ray: 

1 + .^^= 1^"^"|^ (1.24) 



A fluid is modelled as a continuum with a well-defined average 4- velocity field -u", 
where u°'Ua = — 1- This hydrodynamic description requires that the mean collision time 
is much less than any macroscopic characteristic time (such as the expansion time in 
an expanding universe); equivalently, the mean free path must be much less than any 
macroscopic characteristic length. For a perfect fluid, is uniquely defined^ as the 
4-velocity relative to which there is no particle current, i.e. 

= n-u" (1.25) 

where n is the number density. 

The field of comoving observers defines a covariant splitting of spacetime into time 
+ space (1 + 3) via the projection tensor 

hap = gafS + UaUfj = hfSa hafiu'^ = , h^^hf^p = h^p 

/i"a = 3 , Kpq^ = qa if g^w" = (1.26) 

which projects at each point into the instantaneous rest space of the fluid/ observer, and 
provides a 3-metric in the rest space. In the comoving lOF 

= (1, 0) , haf3 = diag (0, 1, 1, 1) , Kpq'^q^ = q- q 

where qaU°' = 0. This allows us to compare relativistic fluid kinematics and dynamics 
with its Newtonian limit. 

The covariant time derivative along is 

A"-p... = A^-p..,X (1-27) 
and describes the rate-of-change relative to comoving observers. In the comoving lOF 

A ,--^A .... 

The covariant spatial derivative is 

D„/ = hjf^p (1.28) 
D,g;3 = ha^hp-'V^q, (1.29) 
D„a/3^ = K^hp^'h^'^V^.a^^, etc. (1.30) 



^If the fluid is out of equilibrium as a result of dissipative effects, then there is no unique average 
4-velocity 
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and describes spatial variations relative to comoving observers. In the comoving lOF, 
with g° = (g°, q) 

D„/ = (0, V/) , D°g„ = V ■ g , e'^^^D^qu = (V x q)' (1.31) 



Any spacetime vector can be covariantly split as 

l^" = v4u" + 5", where A = , 5" = /i"^^^ ^ = (1.32) 

For a rank-2 tensor: 

Vap = AUaUp + BaUp + UcCp + F^p (1.33) 

where A = Vapu'^u^, BaU" = = CaU" and 

Fap = K^hp'V^, ^ F^pu'' = = F^pu^ 

For example, if Va/s = Wa-p, then F^p = DpWa- Now Fap may be further decomposed 
into symmetric and skew parts: 

Fap = F(^aP) + F\^ap] 

In the comoving lOF, the corresponding decomposition of the matrix of components Fij 
is simply 

F={F) + [F] = \(f + F^) + \{F- F^) 
and {F) may be further split into its trace and trace-free parts: 

(F) = {ltr f]i+{F) 

The covariant expression of this is 

F{aP) = haP + F^aP> 

where the symmetric, spatial trace-free part of any tensor is defined by 

V<ap> = K^hp'' {\/(^,) - iK./i'^'' (1.34) 
Thus we can rewrite the decomposition ( |1.33|) in the covariant irreducible form 



Vap = Au^U^ + BaUp + UaCp + IV^uh^^h^p + K;a/3> + V[^,u\h^ ^K" p (1.35) 



Now we are ready to define the quantities that covariantly describe the fluid kinematics. 
These quantities are simply the irreducible parts of the covariant derivative of the fluid 
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4-velocity. With Vaf3 = Ua;f3, we have A = = Ca since Ua^pu"" = 0, and then Ba 
—Ua^pu^ = —Ua- Thus (|1.35|) gives 



Ua;(3 = Hha/B + (Ja/B + ^a/B — UaUj3 where 3H = U°' -a = D"Uq, , 

(^al3 = U<a;l3> = D</3«a> , ^^a/? = hj" h Uy^-y] = D[^Uq] (1.36) 

In a comoving lOF at a point P, v is zero at P, but its derivatives are not, and we find 
using (|1.31|) that 

3H = V-v, e'^^ujjk = - {V X vy 

so that H generahses the Newtonian expansion rate and ujais generahses the Newtonian 
vorticity. Similarly, it can be seen that cTq,/? is the relativistic generalisation of the Newto- 
nian shear. These kinematic quantities therefore have the same physical interpretation as 
in Newtonian fiuids. A small sphere of fiuid defined in the lOF of a comoving observer at 
t = 0, and then measured in the observer's lOF a short time later, undergoes the following 
changes: 

• due to H, its volume changes but not its spherical shape; 

• due to cTq,/3, its volume is unchanged but its shape is distorted in a way defined by 
the eignevectors (principal axes) of the shear; 

• due to uJai3, its volume and shape are unchanged, but it is rotated about the direction 
V X V. 

The expansion rate defines a comoving scale factor a that determines completely the 
volume evolution: 

H = - (1.37) 
a 

1.3 Conservation Laws - Perfect Fluids 

Assuming there are no unbalanced creation/ annihilation processes, particle number is 
conserved in the fiuid. In an lOF, this is expressed via the continuity equation 

dn - , ^, 
_ + V ■ (nv) = 



By ( |1.25|) , the covariant form of particle conservation is 

n'^^a = -v^ h + 3Hn = -v^ na^ = comoving const (1.38) 

where ( |1.37| ) was used to show that the comoving particle number N oc na^ is constant. 

A perfect fiuid is described by its 4-velocity m", number density n, energy (or mass- 
energy) density p, pressure p and specific entropy 5*. In a comoving lOF, the pressure 
is isotropic and given by the Newtonian stress tensor Tij = p6ij. This can be covariantly 
combined with the energy density into the symmetric energy-momentum tensorQ 

Tap = PUaUp + pKp (1.39) 



^The form of the energy-momentum tensor may be justified via relativistic kinetic theory 
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so that Too = P = energy density, Tij = Tij = momentum density, Toi = 0. Just as the 



divergence of produces a conservation law (|1.38| ), so too does the divergence of T°^: 



rpa/B 



;/3 



p + 3H{p + p) = 

{P + P)Ua + DaP = 



;i.4o) 
;i.4i) 



In a comoving lOF these become 



dp 
dt 



+ {p + p)V-v = 0, (p + p) 



dv 
di 



-Vp 



so that ( |1.40|) is an energy conservation equation, generahsing the mass conservation 
equation of Newtonian fluid theory, while ( |1.41| ) is a momentum conservation equation, 
generalising the Euler equation. (In relativity, the pressure contributes to the effective en- 
ergy density.) The energy-momentum conservation equation also follows from Einstein's 
field equations (|1.18|) and the contracted Bianchi identity ( |1.19| ). Equivalently, the con- 
servation equation ensures that the identity holds, i.e. that this integrability condition of 
the field equations is satisfied. 

Finally, the entropy is also conserved. In a comoving lOF, there is no entropy flux, 
and the specific entropy S is constant for each fluid particle. The covariant expression of 
this statement is 

5";a = where = ^ 5 = (1.42) 



where ( |1.38D was used. Note that S is constant along fluid particle worldlines, and not 
throughout the fluid in general. If S is the same constant on each worldline - i.e. if 
= as well as S" = 0, so that S^a = - then the fluid is called isentropic. 



1.4 Equilibrium Thermodynamics 

A perfect fluid is characterised by (n", 5", T"^), or equivalently by (n, p,p, 5, u"), subject 
to the conservation laws above. What are the further relations amongst the thermody- 
namic scalars n, p, p, 5* and T, the temperature? Firstly, the temperature is deflned via 
the Gibbs equation 

TdS = d(f^\+pd(-\ (1.43) 



where df = f^adx"'- Secondly, thermodynamical equations of state are needed in order to 
close the system of equations. Equations of state are dependent on the particular physical 
properties of the fluid, and are deduced from microscopic physics (i.e. kinetic theory and 
statistical mechanics), or from phenomenological arguments. In fact, assuming the metric 
is known (and so leaving aside Einstein's fleld equations), there are 7 equations - i.e. 
(|0|), ([LlOD , (|r41|)i, ([Li^ , ([Lisp - for 8 variables - i.e. n, p,p,Ui, S,T. Thus a single 



scalar equation of state will close the system. 

The Gibbs equation shows that in general two of the thermodynamical scalars are 
needed as independent variables. For example, taking n, p as independent, the remaining 
thermodynamical scalars are p{n, p), S{n, p),T{n, p), and given any one of these, say p = 
pin, p), the others will be determined. Often a barotropic equation of state for the pressure 
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is assumed, i.e. p = p{p). By the Gibbs equation, this imphes S is constant (see below), 
i.e. the fluid is isentropic. 

The adiabatic speed of sound Cg in a fluid is given in general by 



For a perfect fluid, this becomes 



= t 
P 



;i.44) 



1.45) 



as can be seen by choosing p, 5* as independent variables, and using the fact that 5 = 0: 

S 



■ dp\ . 



^ dp\ 



The preceding considerations are phenomenological and mathematical. If the fluid 
model is based on microscopic physics, further conditions are imposed. For example, if 
the fluid is a collision-dominated gas in equilibrium, then relativistic kinetic theory, based 
essentially on imposing energy-momentum conservation at a microscopic level, leads to 
stringent conditionsFI. If m > is the rest mass of the particles and 



m 



m 
T 



then the following conditions hold:| 



Pm = 



mn = Cq- 



KM 

P 



p = nT 



P = Co 



K,{P) ^_^K,iP) 



;i.46) 
;i.47) 

1.48) 



P ' ^ P' 

where Cq is a constant and are modified Bessel functions of the second kind. Further- 



more, 



1^ 1.46 ) shows that is a Killing vector field, so that the spacetime is stationary. 
In particular, (|1.36| ) implies 



H = 0, 

and then ([OD, ([Tig) lead to 



Ur 



-D„ In T , a. 



n=p=p=T=0 



1.49) 



1.50) 



Thus if the perfect fluid is a relativistic Maxwell-Boltzmann gas in equilibrium, severe 
restrictions are imposed not only on the fluid dynamics but also on the spacetime geom- 
etry. 

"^Note that kinetic theory incorporates assumptions about the interactions of particles, in particular 
that the interactions are described by the Boltzmann collision integral. 
^See j^. In standard units, (3 — mc^/kT. 
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In the case of a gas of massless particles in coUisional equihbrium, the conditions are 
less severe: 



T 



n 



p = 3p 



H = 
3nT 



T 

-7^ , (^af3 



;i.5i) 
;i.52) 



Thus is a conformal Killing vector field, so that expansion is possible in equilibrium. 

Kinetic theory shows that a purely phenomenological approach to fluid thermodynam- 
ics holds potential problems in the form of hidden or unknown consistency conditions that 
may be violated. Any phenomenological model needs to be applied with caution. 

The best motivated barotropic perfect fluid model is that for incoherent radiation or 
massless particles, for which p = |p, as in ( |1.52D . The energy conservation equation ( |1.40| ) 
integrates, on using (|1.37|) : 

p = (comoving const)a~^ (1.53) 
Cold, non-relativistic matter is often modelled as pressure-free 'dust', so that 



p 







P 



(comoving const) a 



mn 



;i.54) 



A kinetic theory motivation for the dust model arises from (|1.47| ), (|1.48|) in the limit 

/3> 1: 

p = nT , p ^ mn + |nT where T m (1.55) 

The energy density is p ~ n{mc^ + e), where mc^ is the rest mass energy per particle, 
and e = ^kT is the thermal energy per particle. While (|1.55 ) is still reasonable at high 
temperatures (e.g. for the electron, m ~ lO^K, and ( |1.55[ ) should be very accurate for T 
up to about lO^K), the exact limiting dust case is only reasonable at low temperatures, 
when random velocities are negligible. Of course the hydrodynamic description is no 
longer valid in this limit. 

We can find the evolution of the temperature easily in the case of radiation. Comparing 
( |L5^ ) and ([1331) , we get 

radiation: T oc - (1.56) 

a 

In the general case, the Gibbs equation ( |1.43[ ) can be written as 

'p+p 



dS = - 



Tn^ 



i\ 1 
■]dn+ —dp 
J In 



and the integrability condition 



becomes 



dTdn dndT 

n— + ( +p)^ = T^ 
dn dp dp 



;i.57) 



Furthermore, since T = T{n, p), it follows on using number and energy conservation (|1.38|) 
and ( ITiq ) that 

r dT dT' 
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and then ( |1.57D implies 

^-(11 

From the derivation of ( |1.58|) , we see that it will hold identically if the Gibbs integrability 
condition, number conservation and energy conservation are satisfied. 

This equation holds for any perfect fluid. For non-relativistic matter (|1.55|) gives 

p = |(p — rnn) 

so that (|1.58|) implies: 

non-relativistic matter: T oc — (1.59) 

This shows that the mean particle speed decays like a~^, since the thermal energy per 
particle is e ~ |fcT ^ |mt;^. Strictly, the limiting case of dust has T = 0, but if dust 
is understood as negligible pressure and temperature rather than exactly zero pressure, 
then (iLSQl) holds. 

Note that the Gibbs integrability condition shows explicitly that one cannot inde- 
pendently specify equations of state for the pressure and temperature. This is clearly 
illustrated in the barotropic case. 



Barotropic Perfect Fluids 



With p,p as the independent variables in the Gibbs equation ( |1.43| ) in the general 
perfect fluid case, we flnd: 



(P + P) 



dS 



dn dn 
—dp + —dp 
op op 



+ 



n 



(p + p) 



dp 



n 



dn 



p + p dp 
n 



dp - ^dp 
dp 



h pdn 
p + p p p dp 
pdn ^ 
p dp dp 



dp - ^dp 



where we used the conservation equations (|1.38|) and (|1.40|) . Thus, for any perfect fluid 



n'^TdS = ip + p) 



dn 
dp 



p 



-dp — dp 



P 



;i.60) 



Suppose now that the pressure is barotropic: p = p{p). It follows immediately from ( |1.60| ) 
that dS = 0, i.e. the fluid is isentropic. 

The same conclusion follows in the case of barotropic temperature. Choosing p, T as 
the independent variables, we flnd 



n^TdS =ip + p) 



dn 
df 



T 
— ( 

P 



-dp — dT 
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so that T = T{p) implies dS = 0. 

If the pressure and temperature are barotropic, then the Gibbs integrabihty condition 
( |1.57| ) strongly restricts the form of T(p): 



P 



pip) and T = Tip) =^ T oc exp / J* — (1-61) 

J p{p) + p 



The radiation and dust models are cases of a linear barotropic equation of state that 
is often used for convenience 

p = (7 — l)p =^ p = (comoving const)a~^''' (1-62) 

By ( p. .451) , the speed of sound is Cg = 1/7 — 1. For fluids which have some basis in kinetic 
theory, one can impose the restriction 1 < 7 < |. In principle | < 7 < 2 still leads 
to an allowable speed of sound (7 = 2 is known as 'stiff matter'). The false vacuum of 
inflationary cosmology may be formally described by the case 7 = 0. 
If ( p..62|) holds then the Gibbs integrabihty condition (|1.57| ) becomes 



dT dT 
''^^^P^p=^^-^^^ 

whose solution by the method of characteristics yields 



;i.63) 



where F is an arbitrary function. By ( |1.38| ) and ( p.. 621) , F is a comoving constant, i.e. 



F = 0. If T is also barotropic, then F is constant and we have a power-law form with 
fixed exponent for the temperature: 

T oc p(^-i)/^ (1.64) 



The same result follows directly from ( |1.61| ). 



Note that (|1.62|) and (|1.64| ) are consistent with the ideal gas law p = riT. For dissipa- 



tive fluids, this is no longer true. 



1.5 Example: Cosmological Fluids 

The Ricci identity ( p..l6| ) for the fluid 4-velocity, appropriately projected and contracted, 
together with the field equations ( |1.18|) , leads to an evolution equation for the expansion 
rate 

3if + 3/72 _ u^,^ + a^pa''^ - uj^pu""^ = -l(p + 3p) (1.65) 

known as Raychaudhuri's equation. 

In the standard FRW cosmological models, the rest spaces of comoving observers mesh 
together to form spacelike 3-surfaces {t = const}, where t is proper time for comoving 
observers. Each comoving observer sees that there are no preferred spatial directions - i.e. 
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the cosmic 3-surfaces are spatially isotropic and homogeneous. Thus for any covariant 
scalar / and vector Va 



and 



D„/ = / = f{t)] , h/vp = [^v^ = V{t)ua\ 

Ua;f3 (X haf^ <^ Ma = , (Tap = = 



Raychaudhuri's equation (|1.65|) reduces to 

3H + 3H^ = -|(p + 3p) 



;i.66) 



The momentum conservation equation ( [L.41| ) is identically satisfied. Since p = p{t), 
p = pit), it follows that p = p{p), i.e. one may assume a barotropic equation of state (for 
a single-component fiuid). Then ( |1.66D and the energy conservation equation ( |1.4(j| ) are 
coupled equations in the 2 variables H,p, and can be solved for a given p{p). However 
it is more convenient to use the Friedmann equation, the (0, 0) field equation, which is a 
first integral of the Raychaudhuri equation: 



;i.67) 



where a{t) is the scale factor defined by (|1.37|) and k = 0, ±1 is the curvature index for 
the cosmic 3-surfaces, which by symmetry are spaces of constant curvature. In comoving 
spherical coordinates, the FRW metric and 4-velocity are 



dr"^ 



1 — kr"^ 



u 



f 



:i.68i 



where dVl^ is the metric of the unit sphere. 

The expansion of the universe {H > 0) is confirmed by the systematic redshift in 
electromagnetic radiation that reaches us from distant galaxies. By ( |1.24| ) and (|1.68|) 



1 + z 



a{t 



R) 



a{ti 



;i.69) 



showing that a is increasing, so that by ( |1.62|) p is decreasing. The early universe is very 
hot, as confirmed by the after-glow we observe in the form of the cosmic microwave back- 
ground radiation. The early universe is modelled by a radiation fiuid ( |1.53D , while the late 
universe is cold and the dust model ( |1.54| ) is appropriate. The transition from radiation- 
to matter-domination requires a careful analysis, and has to deal with the interaction 
between radiation and matter. This covers the recombination era of the universe, and 
involves dissipative processes which I will discuss later. 

Leaving aside this transition (which occupies a very short time in the evolution of 
the universe), the matter and radiation are effectively non-interacting. In the super- 
hot conditions of the early universe, matter particles are ultra-relativistic and effectively 
massless, so that a radiation fluid in equilibrium is a good approximation. In the late 
universe, ( 1.53| ) and ( 1.54 ) show that the energy density of radiation is negligible compared 
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to that of matter, and the dust model becomes reasonable. For the flat universe case 
(A; = 0), ([L53|) , ([L54|) and {^^) lead to the solutions: 



radiation: a oc t^^'^ , matter: a oc t^^^ (l-^O) 

Einstein's theory predicts that a radiation FRW universe will begin at t = with 
infinite energy density and temperature. However, for times less than the Planck time 
tp ^ 10~^^ sec, quantum gravity effects are expected to become dominant, and Einstein's 
theory will no longer hold. As yet, no satisfactory quantum gravity theory has been 
developed, and models of the very early universe are necessarily speculative. One fairly 
successful model, which applies during the semi-classical period between the quantum era 
and the classical Einstein era, is inflation. Inflationary models aim to answer some of the 
problems that arise in the standard classical cosmology (the 'big bang' model). 

In these models, the energy density of the universe is dominated by a scalar field at 
around 10~^^ — 10~^^ sec. The pressure of the scalar field is negative, which acts like an 
effective repulsive force, leading to accelerated expansion, or inflation, during which the 
scale factor a increases by around 10^^. Although the scalar field is not a fluid, it has an 



energy-momentum tensor of the perfect fluid form ( |1.39|) . The condition for accelerated 
expansion is a > 0, so that, by 



inflation ^ d > ^ p < -|p (1-71) 
Particular forms of inflation are exponential inflation in a flat FRW universe, for which 

a oc exp{Hit) and p = —p (1-72) 
and power-law inflation, for which a oc t^, > 1. 
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Chapter 2 

Dissipative Relativistic Fluids 



Perfect fluids in equilibrium generate no entropy and no 'frictional' type heating, because 
their dynamics is reversible and without dissipation. For many processes in cosmology 
and astrophysics, a perfect fluid model is adequate. However, real fluids behave irre- 
versibly, and some processes in cosmology and astrophysics cannot be understood except 
as dissipative processes, requiring a relativistic theory of dissipative fluids. 

In order to model such processes, we need non-equilibrium or irreversible thermody- 
namics. Perhaps the most satisfactory approach to irreversible thermodynamics is via 
non-equilibrium kinetic theory. However, this is very complicated, and I will take instead 
a standard phenomenological approach, pointing out how kinetic theory supports many 
of the results. A comprehensive, modern and accessible discussion of irreversible thermo- 
dynamics is given in [Q. This text includes relativistic thermodynamics, but most of the 
theory and applications are non-relativistic. A relativistic, but more advanced, treatment 
may be found in Q (see also 0, 0). 

Standard, or classical, irreversible thermodynamics was first extended from Newtonian 
to relativistic fluids by Eckart in 1940. However, the Eckart theory, and a variation of 
it due to Landau and Lifshitz in the 1950 's, shares with its Newtonian counterpart the 
problem that dissipative perturbations propagate at infinite speeds. This non-causal 
feature is unacceptable in a relativistic theory - and worse still, the equilibrium states in 
the theory are unstable. 

The problem is rooted in the way that non-equilibrium states are described - i.e. 
via the local equilibrium variables alone. Extended irreversible thermodynamics takes its 
name from the fact that the set needed to describe non-equilibrium states is extended 
to include the dissipative variables. This feature leads to causal and stable behaviour 
under a wide range of conditions. A non-relativistic extended theory was developed by 
MuUer in the 1960's, and independently a relativistic version was developed by Israel and 
Stewart in the 1970's. The extended theory is also known as causal thermodynamics, 
second-order thermodynamics (because the entropy includes terms of second order in the 
dissipative variables), and transient thermodynamics (because the theory incorporates 
transient phenomena on the scale of the mean free path/ time, outside the quasi-stationary 
regime of the classical theory). 

In this chapter I will give a simple introduction to these features, leading up to a 
formulation of relativistic causal thermodynamics that can be used for applications in 
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cosmology and astrophysics. 



2.1 Basic Features of Irreversible 
Thermodynamics 

For a dissipative fluid, the particle 4-current will be taken to be of the same form as 
( |1.38| ). This corresponds to choosing an average 4-velocity in which there is no particle 



flux - known as the particle frame. At any event in spacetime, the thermodynamic state 
of the fluid is close to a fictitious equilibrium state at that eventQ The local equilibrium 
scalars are denoted n, p,p, S, T, and the local equilibrium 4-velocity is u^. In the particle 
frame, it is possible to choose u'^ such that the number and energy densities coincide 
with the local equilibrium values, while the pressure in general deviates from the local 
equilibrium pressure: 

n = n, p = p, p = p + Il (2.1) 

where U = p — p is the bulk viscous pressure. From now on I will drop the bar on the 
equilibrium pressure and write p + 11 for the effective non-equilibrium pressure: 

Peff = P + n (P ^ Peff ' P-^P) 

The form of the energy-momentum tensor may be deduced from the equilibrium form 
( |1.3y| ) and the general covariant decomposition ( |1.35| ), given that Tq,^ is symmetric: 



Tal3 = PUaU/3 + {P + n)/l«^ + QaUp + qpU^ + TT^p (2.2) 

where 

In a comoving lOF, g^, = (0,g) and tTq,^ = TCij6aS^^, so that g is an energy flux (due to 
heat flow in the particle frame) relative to the particle frame, while vTjj is the anisotropic 
stress. 

Both the standard and extended theories impose conservation of particle number and 
energy-momentum: 

n°;„ = , T"^;/3 = 



Particle number conservation leads to the same equation (|1.38| ) that holds in the equilib- 
rium case. However the equilibrium energy and momentum conservation equations ( 1.40 ) 
and ( |1.41| ) are changed by the dissipative terms in (|2.2|) : 



p + 3i/(p + p + n) + D"g, + 2M„g° + a,^7r°^ = (2.3) 
{p + p + Il)ua + D«(p + n) + D^TTap + u^T^ap 

+ h/qp + {mKp + a^^ + uj^^)q^ = (2.4) 

In irreversible thermodynamics, the entropy is no longer conserved, but grows, ac- 
cording to the second law of thermodynamics. The rate of entropy production is given by 

^Note that the equihbrium states are different at different events, and therefore not subject to differ- 



ential conditions such as ( i.46 ) - ( f .52 ) 
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the divergence of the entropy 4-current, so that the covariant form of the second law of 
thermodynamics is 

S'^-a > (2.5) 
5*" no longer has the simple form in (|1.42| ), but has a dissipative term: 

= Snu'' + ^ (2.6) 

where S = S and T = T are still related via the Gibbs equation ( p..43| ) .p| 

The dissipative part i?" of is assumed to be an algebraic function (i.e. not con- 
taining derivatives) of and T"'^, that vanishes in equilibrium: 

= R'{nf^^ T'"') and = 

This assumption is part of the hydrodynamical description, in the sense that non-equilibrium 
states are assumed to be adequately specified by the hydrodynamical tensors n°',T°'^ 
alone.0 The standard and extended theories of irreversible thermodynamics differ in the 
form of this function. 



2.2 Standard Irreversible Thermodynamics 

The standard Eckart theory makes the simplest possible assumption about i?" - i.e. that 
it is linear in the dissipative quantities. The only such vector that can be algebraically 
constructed from (11, qa, Ha/s) and is 

/(n,p)nM" + (/(p,n)g° 

Now the entropy density —UaS°' should be a maximum in equilibrium, i.e. 

" d 







eqm 



This implies / = 0. In a comoving lOF, g^/T = (0, q/T), which is the entropy flux due 
to heat flow. Thus g = 1 and ( |2.6|) becomes 



= Snu'' + ^ (2.7) 

Using the Gibbs equation ( |1.43| ) and the conservation equations (|1.38| ) and (|2.3[), the 
divergence of (|2.7|) becomes 



SHU + (D„ In T + M„) q" + a„^7r 



«/3 



(2.^ 



^In extended thermodynamics, this is the Israel-Stewart approach. An alternative approach is to 
extend the Gibbs equation by including dissipative terms, and to use a generalised temperature, specific 
entropy and pressure. The two approaches agree near equilibrium |Q. 

■^In kinetic theory, this corresponds to truncating the non-equilibrium distribution function - via the 
Grad 14-moment approximation H. 
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Notice that the equihbrium conditions ( p.. 491 ) from kinetic theory lead to the vanishing of 



each factor multiplying the dissipative terms on the right, and therefore to S^'-a = 0. 

From ( p.8|) , we see that the simplest way to satisfy ( p.5|) is to impose the following 
linear relationships between the thermodynamic 'fluxes' 11, qa, Tia/s and the corresponding 
thermodynamic 'forces' H, iia + D„ In T, 0"^/?: 

n = -^CH (2.9) 

g„ = -A(D,T + TnJ (2.10) 

T^ap = -'2r]aai3 (2-11) 

These are the constitutive equations for dissipative quantities in the standard Eckart 
theory of relativistic irreversible thermodynamics. They are relativistic generalisations of 
the corresponding Newtonian laws: 

n = -SCV-tT (Stokes) 
q = — AVT (Fourier) 
TVij = —2'qaij (Newton) 

This is confirmed by using a comoving lOF in ( |2.9|) - (|2.11|) - except that in the relativistic 
case, as discovered by Eckart, there is an acceleration term in ( p.lO| ) arising from the 



inertia of heat energy. Effectively, a heat flux will arise from accelerated matter even in 
the absence of a temperature gradient. 

The Newtonian laws allow us to identify the thermodynamic coefficients: 

• ({p,n) is the bulk viscosity 

• X{p,n) is the thermal conductivity 

• ri[p, n) is the shear viscosity 

Given the linear constitutive equations ( p^) - ( p.ll| ), the entropy production rate 



becomes 



which is guaranteed to be non-negative provided that 

C>0, A>0, r/>0 
Note that the Gibbs equation (|1.43|) together with number and energy conservation 



( |1.38| ) and (|2.3|), leads to an evolution equation for the entropy: 

TnS = - - M«g" - (T„/3vr"'^ (2.13) 

Many, probably most, of the applications of irreversible thermodynamics in relativity 
have used this Eckart theory. However the algebraic nature of the Eckart constitutive 
equations leads to severe problems. Qualitatively, it can be seen that if a thermodynamic 
force is suddenly switched off, then the corresponding thermodynamic flux instantaneously 
vanishes. This indicates that a signal propagates through the fluid at infinite speed, 
violating relativistic causality.0 

'''Even in the Newtonian case, infinite signal speeds present a problem, since physically we expect the 
signal speed to be limited by the maximum molecular speed. 
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2.3 Simple Example: Heat Flow 



For a quantitative demonstration, consider the flow of lieat in a non-accelerating, non- 
expanding and vorticity-free fluid in flat spacetime, where the comoving lOF may be 
chosen as a global orthonormal frame. In the non-relativistic regime the fluid energy 
density is given by ( |1.55| ), and then the energy conservation equation ( p.3|) gives 

9T ^ 



since dn/dt = by ( p.. 381 ). The Eckart law ( ^.1(J| ) reduces to 

q = -AVT 

Assuming that A is constant, these two equations lead to 

dT 2A 

— = XV^T where x = ^ (2-14) 

which is the heat conduction equation. This equation is parabolic, corresponding to 
infinite speed of propagation. 

Apart from causality violation, the Eckart theory has in addition the pathology of 
unstable equilibrium states. It can be argued that a dissipative fluid will very rapidly 
tend towards a quasi-stationary state that is adequately described by the Eckart theory. 
However, there are many processes in which non-stationary relaxational effects dominate.^ 
Furthermore, even if the Eckart theory can describe the asymptotic states, it is clearly 
unable to deal with the evolution towards these states, or with the overall dynamics of 
the fluid, in a satisfactory way. 

Qualitatively, one expects that if a thermodynamic force is switched off, the corre- 
sponding thermodynamic flux should die away over a finite time period. Referring to the 
heat flow example above, if VT is set to zero at time t = 0, then instead of q{t) = for 
t > 0, as predicted by the Eckart law, we expect that 

q{t) = go exp (^-^ 

where r is a characteristic relaxational time for transient heat flow effects. Such a relax- 
ational feature would arise if the Eckart-Fourier law were modified as 

r^+g = -AVT (2.15) 

This is the Maxwell-Cattaneo modification of the Fourier law, and it is in fact qualitatively 
what arises in the extended theory. 

With the Maxwell-Cattaneo form ( p.l5|) , the heat conduction equation (p.l4|) is mod- 



ified as 

^For examples and further discussion, see H, [p|- 
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which is a damped wave equation. A thermal plane-wave solution 

T oc exp i{k ■ x — ut) 

leads to the dispersion relation 



X 



so that the phase velocity is 

V 



UJ 



+ iuj 



1/2 



In the high frequency limit, i.e. uj ^ r ^, we see that 

The high-frequency limit gives the speed of thermal pulses - known as second sound - 
and it follows that this speed is finite for r > 0. Thus the introduction of a relaxational 
term removes the problem of infinite propagation speeds. 

The intuitive arguments of this section form an introduction to the development of 
the extended theory of Israel and Stewart. 



2.4 Causal Thermodynamics 

Clearly the Eckart postulate ( p.7|) for i?" is too simple. Kinetic theory indicates that in 
fact R°' is second-order in the dissipative fluxes. The Eckart assumption, by truncating 
at first order, removes the terms that are necessary to provide causality and stability. The 
most general algebraic form for i?" that is at most second-order in the dissipative fluxes 
is 

^a^^a^ (2.17) 

where /3yi(p, ra) > are thermodynamic coefficients for scalar, vector and tensor dissipa- 
tive contributions to the entropy density, and a^(p, n) are thermodynamic viscous/ heat 
coupling coefficients. It follows from (|2.17|) that the effective entropy density (measured 
by comoving observers) is 

- u,S^ = - ^ (/^on^ + (3iq,q^ + /527r^.7r^'^) (2.18) 

independent of ao, a\. (Note that the entropy density is a maximum in equilibrium.) 
For simplicity, I will assume 

ao = = «! i.e. no viscous/ heat coupling (2-19) 
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This assumption is consistent with hnearisation in a perturbed FRW universe, since the 
couphng terms lead to non-hnear deviations from the FRW background. However, the 
assumption ( p.l9| ) may not be reasonable for non-uniform stellar models and other situ- 
ations where the background solution is inhomogeneous. 

The divergence of the extended current ( p.l7|) - with ( p. 19 ) - follows from the Gibbs 
equation and the conservation equations ( |1.38| ), ( |2.3| ) and 



-n 



?,H + /5on + iT 



n 



Da In T + Ma + I3iqa + \T 



T 



— vr 



T 



-u 



IT. 



[2.20) 



The simplest way to satisfy the second law of thermodynamics ( p.5|) , is to impose, as in 
the standard theory, linear relationships between the thermodynamical fluxes and forces 
(extended), leading to the following constitutive or transport equations^: 



tqU + U = -3CH 



f ^0 c 



-A (D„r + Tm«) 



rjT 



^2 , 



n 

I ( ^1 



VAT' 



;/3 



2r]T 



vr, 



q:/3 



where the relaxational times T^(p, n) are given by 



:2.21) 



[2.22\ 



(2.23) 



(2.24) 



With these transport equations, the entropy production rate has the same non-negative 
form ( |2.12| ) as in the standard theory. 

Because of the simplifying assumption (|2.19|) , there are no couplings of scalar/ vector/ 
tensor dissipative fluxes. As well as these viscous/ heat couplings, kinetic theory shows 
that in general there will also be couplings of heat flux and anisotropic pressure to the 
vorticity - which, unlike the shear, does not vanish in general in equilibrium (see ( |1.49| )). 
These couphngs give rise to the following additions to the right hand sides of ( |2.22| ) and 
( p.23| ) respectively: 

+AT7iu;Q/3g^ and + 2r]^2TT^ <a^fi>iL 

where 7i(p, n), 72(p, n) are the thermodynamic coupling coefficients. In a comoving lOF, 
( |1.31| ) shows that the addition to ( ^^.22] ) has the form 



AT711I; X q where 



^This linear assumption is in fact justified by kinetic theory, which leads to the same form of the 
transport equations 
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If the background solution has zero vorticity, as is the case in a perturbed FRW universe, 
then these vorticity coupling terms will vanish in linearised theory. However, they would 
be important in rotating stellar models, where the background equilibrium solution has 

The terms in square brackets on the right of equations ( p. 21 ) - ( 2.23|) are often omit- 



ted. This amounts to the implicit assumption that these terms are negligible compared 
with the other terms in the equations. I will call the simplified equations the truncated 
Israel-Stewart equations. One needs to investigate carefully the conditions under which 
the truncated equations are reasonable. This will be further discussed in the next chap- 
ter. The truncated equations, together with the no-coupling assumption (|2.19|) , are of 
covariant relativistic Maxwell-Cattaneo form: 

Toll + n = -3C^r (2.25) 
Tihjqp + q^ = -A(D„T + Tu„) (2.26) 
TiK^hp'Ti^^ + Tiap = -2r]aaf3 (2.27) 

The crucial difference between the standard Eckart and the extended Israel-Stewart 
transport equations is that the latter are differential evolution equations, while the former 
are algebraic relations. As we saw in the previous section, the evolution terms, with 
the relaxational time coefficients r^, are needed for causality - as well as for modelling 
high-frequency or transient phenomena, where 'fast' variables and relaxation effects are 
important. The price paid for the improvements that the extended causal thermodynamics 
brings is that new thermodynamic coefficients are introduced. However, as is the case 
with the coefficients (, X,7] that occur also in standard theory, these new coefficients may 
be evaluated or at least estimated via kinetic theory. The relaxation times ta involve 
complicated collision integrals. In fact, they are usually estimated as mean collision 
times, of the form 

1 

(2.28) 



nav 

where cr is a collision cross section and v the mean particle speed. 

It is important to remember that the derivation of the causal transport equations is 
based on the assumption that the fluid is close to equilibrium. Thus the dissipative fluxes 
are small: 

|n|<p, (7r„^7r°'3)'^'<p, (g„g")'/'<p (2.29) 

Consider the evolution of entropy in the Israel-Stewart theory. The equation (|2.13|) 
still holds in the extended case: 

TnS = -3HU - - - a^/jvr"^ (2.30) 

Consider a comoving volume of fluid, initially of size ajj, where a is the scale factor defined 
in general by ( |1.37|) . The entropy in this comoving volume is given by 

E = a^nS (2.31) 

Then, by virtue of number conservation ( |1.38| ) and ( p.30| ), it follows that the growth in 
comoving entropy over a proper time interval to ^ t is 

S(t) = So - ^* ^ + q'',^ + ii^q'^ + a^^^vr"^) dt (2.32) 



23 



The second law, which is built into the theory, guarantees that S(t) > Sq. However, it is 
possible that the local equilibrium specific entropy S is not increasing at all times - but 
the effective, non-equilibrium specific entropy —UaS'^/n is monotonically increasing [Q. 

Next we look at the temperature behaviour in causal thermodynamics. The Gibbs 
integrability condition (|1.57 ) still holds: 



dn dp dp 



(2.33) 



However, the change in the energy conservation equation (|2.3| ) leads to a generalisation 



of the temperature evolution ( |1.58| ): 



T 
f 



-3H 



SHU + + + (t^.tt"^ 



(2.34) 



-3H ( ^ 1 +nS 

Note that if the Gibbs integrability condition, number conservation and energy conser- 



vation are satisfied, then the evolution equation ( |2.34| ) will be an identity. This evolution 
equation shows quantitatively how the relation of temperature to expansion is affected by 
dissipation. The first term on the right of ( |2.34| ) represents the cooling due to expansion. 
In the second, dissipative term, viscosity in general contributes to heating effects, while 
the contribution of heat flow depends on whether heat is being transported into or out of 
a comoving volume. 

If instead of {n, p) we choose (n, T) as independent variables, then the Gibbs integra- 
bility condition ( p.33| ) becomes 



rpdp ^ dp 
dT dn 



P + P 



(2.35) 



and the temperature evolution equation ( p.34| ) becomes 

1 



T 
T 



a/3 



(2.36) 



nS 



1 



{dp/dT)r 



Finally, we consider briefly the issue of equations of state for the pressure and tem- 
perature in dissipative fluids. Using the energy conservation equation (^.31), the Gibbs 
equation in the form ( |1.60| ) generalises to 



n'TdS 



nV 



3H{p + p)+V 



-dp — dp 
P 



where the dissipative term is 



V = SHU + + + a^^Tr"^ 



(2.37) 
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It follows that in the presence of dissipation, barotropic pressure no longer forces dS to 
vanish: 



dS=^ 
nT 



V 



dp 



(2.38) 



_3H{p + p)+V_ 

As in the equilibrium case, it remains true, via the Gibbs integrability condition, that 
barotropic T = T{p) together with p = (7 — l)p leads to the power-law form ( 1.64|) 
for the temperature. However, in the dissipative case, these relations are not in general 
compatible with the ideal gas law p = nT: 



p = nT , ^=(7 — l)p, T (X p 



(7-l)/7 



I/t 

n (X p ' ' 
h 1 p 
n 7 p 



Then number and energy conservation imply V = 0. 

However, it is possible to impose the 7-law and the ideal gas law simultaneously, 
provided the temperature is not barotropic. The temperature evolution equation ( |2.36|) 
and energy conservation (|2.3|) give 



T 
T 



7-1 

. 7 , 



p 7p 



V 

nT, 



(2.39) 



These results have interesting implications for a dissipative fluid which is close to a 
thermalised radiation fluid, i.e. p = ^p. If we insist that p = nT, then the Stefan- 
Boltzmann law p (x T^ cannot hold out of equilibrium. Alternatively, if we impose the 
Stefan-Boltzmann law, then the ideal gas law cannot hold unless the fluid returns to 
equilibrium. 
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Chapter 3 



Applications 
Astrophysics 



to Cosmology and 



The evolution of the universe contains a sequence of important dissipativc processes, 
including: 

• GUT (Grand Unified Theory) phase transition (t a; 10~^"^ sec, T !^ 10^^ K), when 
gauge bosons acquire mass (spontaneous symmetry breaking). 

• Reheating of the universe at the end of inflation (at about 10~^^ sec), when the 
scalar field decays into particles. 

• Decoupling of neutrinos from the cosmic plasma (t ~ 1 sec, T ^ 10^° K), when 
the temperature falls below the threshold for interactions that keep the neutrinos in 
thermal contact. The growing neutrino mean free path leads to heat and momentum 
transport by neutrinos and thus damping of perturbations. Shortly after decoupling, 
electrons and positrons annihilate, heating up the photons in a non-equilibrium 
process. 

• Nucleosynthesis (formation of light nuclei) {t ^ 100 sec). 

• Decoupling of photons from matter during the recombination era {t ~ 10^^ sec, T ^ 
10^ K), when electrons combine with protons and so no longer scatter the photons. 
The growing photon mean free path leads to heat and momentum transport and 
thus damping. 

Some astrophysical dissipative processes are: 

• Gravitational collapse of local inhomogeneities to form galactic structure, when 
viscosity and heating lead to dissipation. 

• Collapse of a radiating star to a neutron star or black hole, when neutrino emission 
is responsible for dissipative heat flow and viscosity. 

• Accretion of matter around a neutron star or black hole. 
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Further discussion of such processes can be found in [Q, [|TU] (but not from a causal 
thermodynamics standpoint). The apphcation of causal thermodynamics to cosmology 
and astrophysics remains relatively undeveloped - partly because of the complexity of the 
transport equations, partly because all of the important dissipative processes have been 
throughly analysed using the standard theory or kinetic theory or numerical methods. 

Causal bulk viscosity in cosmology has been fairly comprehensively investigated - 
see illl] - [0. Shear viscosity in anisotropic cosmologies has been considered in 111 



while heat flow in inhomogeneous cosmologies has been discussed in Causal 
dissipation in astrophysics has been investigated in |§ , - ||2^ . In all of these papers, 
it is found that causal thermodynamic effects can have a significant impact and can lead 
to predictions very different from those in the standard Eckart theory. 

In this chapter I will briefly discuss some overall features of causal thermodynamics in 
a cosmological/ astrophysical setting, and then conclude with a more detailed discussion 
of bulk viscosity in an FRW universe, which is the most accessible problem. 



3.1 General Features of Cosmic Dissipation 

The expanding universe defines a natural time-scale - the expansion time = a/ a. 
Any particle species will remain in thermal equilibrium with the cosmic fluid so long as 
the interaction rate is high enough to allow rapid adjustment to the falling temperature. 
If the mean interaction time is tc? then a necessary condition for maintaining thermal 
equilibrium is 

< (3.1) 

Now tc is determined by 

tc = ^ (3.2) 

nav 

where a is the interaction cross-section, n is the number density of the target particles 
with which the given species is interacting, and v is the mean relative speed of interacting 
particles. 

As an example, consider neutrinos in the early universe. At high enough temperatures, 
the neutrinos are kept in thermal equilibrium with photons and electrons via interactions 
with electrons that are governed by the weak interaction. The cross-section is 

(Jn, = goT^ (3.3) 

where Qq is a constant. The number density of electrons is oc T^, by ( |1.52| ), since the 
electrons are effectively massless at these very high temperatures. Since v = 1, ( |3.2|) gives 
tc oc T~^. By ( |1.56| ), we can see that H oc T^. Thus 

tcH=[^)' (3.4) 

and using ( ^.1| ) and the numerical values of the various constants, it follows that the 
neutrinos will decouple for 

T < ^ 10^° K 
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Other cosmic decoupling processes may be analysed by a similar approach. The differ- 
ences arise from the particular forms of cr(T), n(T) and H{T). For example, in the case of 
photons interacting with electrons via Thompson scattering, the Thompson cross-section 
is constant, while the number density of free electrons is given by a complicated equation 
(the Saha equation), which takes account of the process of recombination. The expansion 
rate H is also fairly complicated, since the universe is no longer radiation-dominated. 
One finds that the decoupling temperature is about 10^ K. 

In the case of a collapsing star, similar arguments are applied - except that the char- 
acteristic time in this case is determined by the rate of collapse, which is governed by 
stellar dynamics. For example, for neutrinos in the core of a neutron star, interactions 
with electrons and nucleons determine an interaction time that must be compared with 
the collapse time to estimate the decoupling conditions for the neutrinos - after which 
they transport heat and momentum away from the core. 

The entropy generated in a dissipative process that begins at to and ends at to + 
is given by ( |2.32| ): 



AS = - 



-(3/7n + g";, + n,g" + a,^7r"'3)rft (3.5) 



For example. At could be the time taken for a decoupling process in the universe or a 
star. 

The observed universe has a high entropy, as indicated by the high number of photons 
per baryon, about 10^. This gives a total entropy in the observable universe of about 10®^. 
Inflationary cosmology predicts that nearly all of this entropy is generated by the reheating 
process at the end of inflation - i.e. that all other dissipative processes in the evolution 
of the universe make a negligible contribution to entropy production by comparison. In 
this model, the formula ( |3.5| ) would have to be modified to include the dissipation not 
just from the fluid effects that we have been discussing, but also from particle production. 
Particle production, at a rate i/, leads to non-conservation of particle number, so that 
( |1.38| ) is replaced by 

n°';a = n + 3Hn = un (3.6) 

Then it is found that u contributes to entropy production. The contribution from particle 
production may be modelled as an effective bulk viscosity. 

Many dissipative processes are well described by a radiative fluid - i.e. a fluid con- 
sisting of interacting massless and massive particles. The radiative fluid is dissipative, 
and kinetic theory or fluctuation theory arguments may be used to derive the dissipative 
coefficients in terms of the relaxation times ta_ (which are usually assumed equal to the 
appropriate interaction time tc). The results are collected in the table below. The table 
also includes the case of a relativistic Maxwell-Boltzmann gas - i.e. a dilute monatomic 
gas with high collision rate - in both the ultra-relativistic and non-relativistic limits. The 
local equilibrium energy density and pressure are given by the equations of state ( |1.47| ), 
( |1.48| ) (but not subject to the global equilibrium conditions ( |1.49|) , (|1.50|) ). 
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c 


A 


7/ 


X OjV-J-IOj U 1 V J. J. LLlV-l 

(massless/ massive) 


4roT4r2ro 






Maxwell— Boltzmann gasi 
ultra-relativistic (/5 <^ 1) 








Maxwell-Boltzmann gas: 
non-relativistic (/?:»!) 


i/5-Vo 




pr2 



In the table, /3 is given in standard units by 



1^ 



where m is the mass of the matter particles (usually electrons); tq is the radiation constant 
for photons, and | times the radiation constant for massless neutrinos; T is effectively the 
deviation of p/ p from its pure-radiation value: 

(dp/dT),, 




' {dp/dT)n 

where p, p refer to the pressure and energy density of the radiation/ matter mixture as a 
whole. For example, when the matter is non-relativistic, ( |1.52| ) and ( |1.55| ) show that in 
standard units 

p ^ nkT + |roT^ , p ^ mc^n + |nA;T + r^T^ (3.8) 

where n is the number density of matter. 

Note that for both the radiative fluid and the Maxwell-Boltzmann gas, the bulk viscos- 
ity tends to zero in the ultra-relativistic and non-relativistic limits. Bulk viscous effects 
are greatest in the mildly relativistic intermediate regime, /? ~ 1. This discussed further 
in the next section. 

The radiative fluid and Maxwell-Boltzmann gas are perhaps the best motivated dissi- 
pative fluid models. However, their equations of state and thermodynamic coefficients are 
very complicated, and for the purposes of analytical rather than numerical investigations, 
simplified equations are often assumed. These are usually barotropic: 

p = p{p), T = T{p), C = C(P), A = A(p), V = V{P), rA = TA{p) (3.9) 
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However these assumptions are subject to consistency conditions (as shown earher in 
the case of p and T), and may correspond to unphysical behaviour. Whenever such 
assumptions are made in a model, the consequences should be carefully checked. An 
example is given in the next section. 



3.2 Causal Bulk Viscosity in Cosmology 



I will use the simplest case of scalar dissipation due to bulk viscosity in order to illustrate 
some of the issues that arise in modelling cosmological dissipation via Israel-Stewart 
theory. Furthermore, this case covers the standard cosmological models. If one assumes 
that the universe is exactly isotropic and homogeneous - i.e. an FRW universe ( |1.68| ) - 
then the symmetries show that only scalar dissipation is possible - i.e. Qq = = Hap- In 
this event, the no-coupling assumption ( p.l9|) is automatically fulfilled. 

Bulk viscosity arises typically in mixtures - either of different species, as in a radiative 
fluid, or of the same species but with different energies, as in a Maxwell-Boltzmann gas. 
Physically, we can think of bulk viscosity as the internal 'friction' that sets in due to the 
different cooling rates in the expanding mixture. The dissipation due to bulk viscosity 
converts kinetic energy of the particles into heat, and thus we expect it to reduce the 
effective pressure in an expanding fluid - i.e. we expect 11 < for H > 0. This is 
consistent with 5* > by ( |2.13| ): 

TnS = -3Hn (3.10) 

Any dissipation in an exact FRW universe is scalar, and therefore may be modelled 
as a bulk viscosity within a thermodynamical approach. As I have argued in the previous 
chapter, the Israel-Stewart thermodynamics is causal and stable under a wide range of 
conditions, unlike the standard Eckart theory. Therefore, in order to obtain the best 
thermo-hydrodynamic model with the available physical theories, one should use the 
causal Israel-Stewart theory of bulk viscosity. 

Writing out the full Israel-Stewart transport equation (p.21| ) (using r = tq), we get 



rri + n = -3Ci^ - 



3H 



c 
c 



T 
T 



(3.11) 



n 



A natural question is - what are the conditions under which the truncated form (|2.25|) 
is a reasonable approximation of the full Israel-Stewart transport equation? It follows 
from (|3ll|) that if 

T 

holds, then the additional terms in (|3.11|) are negligible in comparison with 3(H. The 
condition (|3.12| ) is clearly very sensitive to the particular forms of the functions p{n,p), 
({n, p) and r(n, p). The temperature is determined on the basis of these particular forms 



< 1 



(3.12) 



by the Gibbs integrability condition (|2.33|) and the evolution equation ( p.34| )[']: 



T 
f 



- = -3H 




n fdT\ 



(3.13) 



^or equivalently by ( 2.35 ) and ( 2.36 ) 



30 



The second term on the right shows that bulk stress tends to counteract the coohng due 
to expansion. 

For simphcity, suppose that the pressure and temperature are barotropic, with p hnear: 



P = (7 - 1)P 



(3.14) 



This pressure equation is not unreasonable if the local equilibrium state is radiation or 
cold matter. Since the temperature is also barotropic, it then follows from the Gibbs 
integrability condition (|2.33|) that as in the perfect fluid case, T must have the power-law 
form ([T^ ): 

2-o^p(7-i)/7 (3_15) 

Thus there is no freedom to choose the form of T(p) - it is a power-law, with index fixed 
by 7. With these forms of p{p) and T(p), we can see that the temperature evolution 
equation ( p.l3| ) is identically satisfied by virtue of the energy conservation equation ( |2.3| ): 



p + 3H{p + p + U) =0 



(3.16) 



A simple relation between r and ( is found as follows. It is shown in the appendix to 
this chapter that 

— ^ = 4 (3.17) 

where Cb is the speed of bulk viscous perturbations - i.e. the non-adiabatic contribution 
to the speed of sound f in a dissipative fluid without heat flux or shear viscosity. The 
dissipative speed of sound is given by 



ct + cl<l 



(3.18) 



where Cg is the adiabatic contribution ( |1.44|) , and the limit ensures causality. When (|3.14|) 
holds, = 7 — 1, so that 

C6<2-7 

We will assume that Cb is constant, like Cg. 

Putting together the thermodynamic relationships ( |3.14| ), ( |3.15| ) and ( p.l7| ), the full 
transport equation (|3.11|) becomes 



1 + 



n 



7Cb VP, 

where the effective relaxation time and bulk viscosity are 

c 



(3.19) 



1 + 3-fTH ' 1 + 3-fTH 

Now the near-equilibrium condition ( p.29|) with (|3.14|) implies 

|n| < p 



(3.20) 
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and shows that the second term in square brackets in ( ^.191 ) is neghgible. Thus the full 
equation leads to a truncated equation with reduced relaxation time and reduced bulk 
viscosity: 

r.n + n = -3C*H (3.21) 

The amount of reduction depends on the size of r relative to H. If r is of the order of the 
mean interaction time, then the hydrodynamical description requires tH < 1. If tH ^ 1, 
then T^: ^ T and C* ~ C- But if tH is close to 1, the reduction could be significant. 

Although this reduction is based on the simplified thermodynamical relations assumed 
above, it indicates that the validity of the truncated Israel-Stewart equation can impose 
significant conditions. More realistic thermodynamical relations will require numerical 
calculations. In the case of a Maxwell-Boltzmann gas, such calculations show that the 
behaviour of the truncated and full theories can be very different. The conclusion seems 
to be that the full theory should be used, unless one is able to derive explicitly - and 
satisfy - the conditions under which the truncated version is adequate. 

Assuming that the FRW universe is flat, the Friedmann equation (|1.67|) is 

p = (3.22) 

By (|16D and (pl) we get 



n = -2H - S'yH^ (3.23) 
and together with (|3.21|) and ( |3.20D , this leads to the evolution equation for H: 



if + (67 + N)HH + §7 [3(7 - cl) + n]h'^ = (3.24) 

where 

N = {tHY^ (3.25) 

is of the order of the number of interactions in an expansion time. Intuitively, when 
^ 1, the fluid is almost perfect, while when A^ is close to 1, the dissipative effects are 
significant. This is confirmed by ( p.24| ). For A^ ^ 1, the equation reduces to 



H + |7if2 ~ 
with the well-known perfect fluid solution: 

H ' 



37(t - to) 



On the other hand, for A^ close to 1, the second derivative in ( |3.24| ) cannot be neglected, 
and the solutions will show a range of behaviour very different from the perfect fluid - 
and the standard Eckart - solutions. (Note that the Eckart limit r ^ is Cf, 00 by 
( p.lTp ; the causality condition (|3.18|) does not hold.) 



Of course, a complete model requires the specification of A^. Consider the ultra- 
relativistic fluid of the early universe, with a particle species whose growing mean free 
path is giving rise to dissipation, such as the neutrino. Suppose that r ^ t^, where tc is 
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the mean interaction time, and that the interaction cross-section is proportional to 



hke the neutrino's. Then by ( |3.4| ) we get 



iv=ir' 



T 



HJ 



(3.26) 



For T ^ T*, we have N ^ 1, and dissipation is neghgible. But for T close to T^, 
dissipation effects become significant. The evolution equation ( p.24| ) becomes 



H + 



HH + 2 



4 - 3c^ + ( — 



H^ = {} 



(3.27) 



One could try to solve this equation perturbatively, by the ansatz 

1 



H 



+ eHi + 0{e' 



I will briefly discuss the question of bulk viscous inflation. Suppose dissipation in the 
cosmic fluid produced sufficiently large bulk viscous stress to drive the effective pressure 
negative and thus initiate inflationary expansion. By ( |1.71| ), using the effective pressure, 
the condition for inflationary expansion is 

-n>j9+|p (3.28) 

For a fluid, this violates the near-equilibrium condition 

|n| <p 

Thus viscous fluid inflation, if it were physically possible, would involve non-linear ther- 
modynamics, far from equilibrium. The Israel-Stewart theory, as well as other versions 
of extended thermodynamics and also Eckart's standard thermodynamics, are all based 
on near-equilibrium conditions, and cannot be applied to inflationary expansion - unless 
one makes the drastic assumption that the linear theory applies in the strongly non-linear 
regime. 

Furthermore, there are serious physical problems with hydrodynamic inflation (with- 
out particle production^). The point is that under conditions of super-rapid expansion 
- i.e. very small expansion time - the hydrodynamic regime requires even smaller inter- 
action time. It is hard to see how the fluid interaction rate could increase to stay above 
the expansion rate under conditions where fluid particles are expanding apart from each 
other extremely rapidly. 

For a satisfactory model of bulk viscous inflation, one needs: (a) a non-linear gener- 
alisation of the Israel-Stewart transport equation (|3.11| ) ;P| (b) a consistent model of fluid 
behaviour under super-rapid expansion and strongly non-linear conditions. 



^See Q, Q for particle production models. 
•^One possible generalisation is developed in p£ 
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On the other hand, the reheating period at the end of inflation can be modelled by 
near-equilibrium theory, and the expansion rate is no longer inflationary. However, a 
thermodynamic model needs to incorporate particle production.^ 



Finally, for those who like to analyse and solve differential equations]^ more than they 
like physical analysis, I will give the evolution equation of H with mathematically more 
general (but physically no more satisfactory) thermodynamic equations of state. Suppose 
p and T are given as above, by ( |3.14| ) and ( |3.15D , but instead of the relation ( |3.17| ), with 
constant Ct, linking r and C, we assume the barotropic forms 



r oc 



(3.29) 



where r and q are constants. 

Then with (|3.29|) , the (non-truncated) evolution equation (|3.11|) becomes 



if + 3 



)7l HH + aiH-'^m + (g - r - 1 + 7"^) R-^H'^ 



(3.30) 



where ai and a2 are constants. One can find special exact solutions, including exponential 
and power-law inflation, and perform a qualitative dynamical analysis of ( |3.3CI| ), or of sim- 
ilar equations arising from different forms for the equations of state and thermodynamic 
coefficients. 



■^See m. 

^See also I?) - [|o|. 
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3.3 Appendix: Bulk Viscous Perturbations 



A comprehensive analysis of the causality and stability properties of the full Israel-Stewart 
theory has been performed by Hiscock and Lindblom 0. They consider general pertur- 
bations - i.e. n, gQ,,7rQ/3 all nonzero - about a (global) equilibrium in flat spacetime, but 
the results are valid in cosmology for short wavelength perturbations. In this appendix, I 
will extract from their complicated general results the special case of scalar perturbations 
(only n 7^ 0), when remarkably simple expressions can be obtained. 

The characteristic velocities for general dissipative perturbations are given by equa- 
tions (110) - (128) in p. The purely bulk viscous case is 



ao = = «! 



1 1 

A' J2 



0; (3^ = 



C 



(3.31) 



(See ( ^TtD and (|3l|) 
Equation (127) of 



(PI).) 

gives the speed of the propagating transverse modes: 



2 (p + p)a'^ + 2ai + /3i 



2/52[/5i(p + p) 







on using ( p.31|) . This is as expected for scalar sound-wave perturbations. Equation (128) 
governing the speed v = vl oi propagating longitudinal modes becomes, on dividing by 
/3o/32 and setting ao = = ai: 



[A(p + p)-l] + 



2n (dT\ (p + p) fdT\ ^ 



' dp\ 1 
dp 



+ 



-{-\ 

Dividing by jSi and taking the limit jSi — ; 



dp 

00, this gives 
1 

5 (P + p)Po 



n 
T 



' s 

'or 

dn 



/5oJ 
1 2 



' dp\ 
Bp 



(3.32) 



The first term on the right is the adiabatic contribution to f^, and the second term is 
the dissipative contribution c^, as in ( |3.17| ). 

It is also shown in (pp 478-480) that causality and stability require 



n,{\)^{p + p){i-\' 



' dp^ 



+ 



1 



s 



> 



for all A such that < A < 1. This condition is shown to hold for all A if it holds for 
A = 1, leading to the requirement 



c 



{p + p)t 



< 1 



(3.33) 



i.e. < 1, as expected. This establishes ( |3.18|) . 
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These results refine and correct the widely-quoted statement in \T^ that C/pr = 1 is 
required by causality. 
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